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Abstract. 

The measured spectral energy distribution and variability 
time scale are used to determine the radiation and magnetic- 
field energy densities in the relativistic plasma that forms the 
gamma-ray emitting jet in the blazar 3C 279. Assuming that 
protons are accelerated as efficiently as electrons to a max- 
imum energy determined by the size and magnetic field of 
the emitting region, we calculate the emissivity of neutri- 
nos produced by protons that interact with the external radi- 
ation field intercepted by the jet. The external radiation field 
provides the most important target photons for photomeson 
production of high-energy neutrinos in flat spectrum radio 
quasars (FSRQs). Because of photomeson interactions with 
this field, km^ neutrino telescopes are predicted to detect 
>, 0.1-1 neutrinos per year from blazars such as 3C 279. 
BL Lac objects are weaker neutrino sources if, as widely 
thought, their 7-ray emission is due to Compton-scattered 
synchrotron (SSC) radiation. 



1 Introduction 

Blazars are among the most powerful accelerators of rela- 
tivistic particles in nature, as shown by EGRET observa- 
tions of 100 MeV - GeV emission from over 60 FSRQs and 
BL Lac objects. Consequently, blazars represent a potential 
source of high-energy neutrinos (see Gaisser et al. 1995 and 
Halzen 2001 for reviews). The bright, highly variable 7 ra- 
diation from blazars argues for the existence of emission re- 
gions with high photon energy densities and intense nonther- 
mal particle populations, which are needed for efficient pho- 
tomeson production (e.g., Mannheim and Biermann 1992; 
Mannheim et al. 2001). Previous papers have considered 
internal synchrotron photons as targets for high-energy pro- 
ton interactions (e.g., Mannheim 1993). Here we show that 
UV photons from external radiation fields, which have earlier 
been proposed as a target photon source for nonthermal elec- 
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trons to produce Compton-scattered gamma-rays (e.g., Der- 
mer and Schlickeiser 1993; Sikora et al. 1994), also provide 
the most important photon source for photomeson production 
of > 30 TeV neutrinos in FSRQs. 

We use the measured spectral fluxes and variability time 
scales of synchrotron and Compton radiation from the blazar 
3C 279 to determine physical parameters of its gamma-ray 
emission region. Lower limits to the Doppler factor 5 are de- 
fined by the condition that the emitting region be transparent 
to gamma rays. Equipartition arguments and observed fluxes 
of the synchrotron and Compton components are used to de- 
termine the magnetic field strength B in the plasma blob. 
Assuming efficient acceleration of power-law energetic pro- 
tons to a maximum energy limited by the magnetic field and 
size scale of the region, we calculate spectral emissivities of 
pion decay neutrinos formed through photomeson produc- 
tion. Because the properties of the external radiation field 
are independent of the bulk Lorentz factor F of the radiat- 
ing region, neutrino emissivity calculations using duty cycle 
estimates give accurate neutrino flux estimates during low 
7-ray states. We predict detectable fluxes of neutrinos from 
7-ray loud FSRQs using km^ neutrino telescope arrays such 
as IceCube. 



2 Constraints on Source Parameters 

2. 1 Synchrotron and external photon energy densities 

The comoving synchrotron photon energy density is given 
by u'g = Ls/{'2,TTr1c5'^), where Lg is the observed bolomet- 
ric synchrotron luminosity, r^ is the comoving radius of the 
blob, here assumed spherical, and primes denote comoving 
quantities. The ratio of the SSC and synchrotron luminosi- 
ties Lssc/Ls ~ u'g/uB, where ub = B^/Stt. For simplic- 
ity, we consider the regime where the quasi-isotropic scat- 
tered external radiation field dominates the direct disk radia- 
tion field (Dermer and Schlickeiser 1994). 

If the 7-rays are due to external Compton-scattered radia- 



tion, then Lec/Ls — S'^u^xt/uB — {S/r)'^u'^^^/uB (Sikora 
1997; Dermer, Sturner, and Schlickeiser 1997). Combining 
these two expressions gives u'^^^/u'g^ (r/(5)^ [Lc/Lssc) 
= a(r/(5)^ = a when viewing within the beam of the jet 
{6 < r^^). Consequently u'^^^ ^ aLg /(27rr^c(5''), where 
a = Lec/ Lssc The energy of the external photons in the 
comoving frame is e^^j = Ttext- Models taking into ac- 
count SSC and external Compton (EC) components show 
that a complete spectral fit requires synchrotron, SSC and 
EC components (Bottcher 1999; Hartman et al. 2001). Thus 
a -- 0.1-1 for BL Lac objects and a - 1-10 for FSRQs. 

The spectral energy density Uph(e') = nieC^e' n'^{e') of 
the photons in the blob frame is given in terms of the lyF^, 
flux /(e) through the relation 
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where we relate the measured variability time scale t^ar to 
the blob radius r^ through the expression r^ ~ ct^arS /(I + 
z), i«ar(d) is the observed variability time scale in days, and 
z is the redshift. The expression e' = e(l + z)/5 relates the 
photon energies measured in the source and observer frames. 

2.2 Magnetic field estimates 

The magnetic field in the blob can be estimated by noting 
that the total soft photon energy density u' ^ = u'^ + u'^^^ ~ 
Ub{Lc/Ls), where Lq ~ Lec+Lssc is the total Compton 
power. Hence ub = (Ls/Lc)u'p^ = {Ls/Lc){l + a)u'^. 
Using the earlier expressions for u'^ and rt gives 
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An independent estimate for B can also be made by intro- 
ducing the equipartition parameter -q — u'^^/u-q for the ratio 
of nonthermal electron to magnetic energy densities in the 
jet. Thus rj ^ 1 corresponds to equipartition, though with 
electrons only. For the measured synchrotron flux density 
Si, ex j/^"' with a^ ~ 0.5 siv ^ 10^^ Hz, this gives 
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in Gauss, where 5'( Jy) is the flux density measured at 10^^i/i2 
Hz. 

2.3 77 transparency constraints on Doppler factor 

Transparency to 77 pair-production attenuation requires that 

r^^(e') = 2(jTn',j(2/e')rfc/(3e') < 1, which can be rewrit- 
ten as a lower limit on 5 using equation (1). 

2.4 Photopion and neutrino production 

Energy losses of relativistic protons (and neutrons) are calcu- 
lated on the basis of standard expressions (see, e.g., Berezin- 
skii and Grigoreva 1988) for the cooling time of relativistic 



protons due to photopion production in p7 collisions. If the 
ambient photons have spectral density n' ^^{e'), then 
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where 7^ is the proton Lorentz factor, e^ is the photon energy 
in the proton rest frame, cr(er) is the photopion production 
cross-section, eth ~ 150 McV is the threshold energy for 
the parent photon in the proton rest frame, and K-p^{e^) is 
the inelasticity of the interaction. The latter increases from 
Ki « 0.2 at energies not very far above threshold where the 
single pion production channels dominate, to K2 ^ 0.5 — 0.6 
at larger values of e^ where multi-pion production dominates 
(e.g. see Berezinskii and Grigoreva 1988, Miicke et al. 1999). 

A detailed recent study of this photohadronic process has 
been given by Miicke et al. (1999). In order to simplify cal- 
culations, we approximate the cross-section cr(er) as a sum of 
2 step-functions CTi(er) and (T2(er) for the (total) single-pion 
and multi-pion channels, respectively, with cti = 380 /ib for 
200 MeV < Er < 500 MeV and cti = outside this region, 
whereas (J2 — 120 /ib at e^ > 500 McV. The inelasticity 
is approximated as K-p^ — Ki and Kp~^ — K2 below and 
above 500 MeV. The results of our calculations in this simpli- 
fied approach are found to give good agreement with earlier 
calculations (Berezinskii and Grigoreva 1988). In particular, 
the time scales for photopion interactions of ultra-relativistic 
cosmic rays with the CMBR are accurately reproduced (see 
also Stanev et al. 2000). This approach also works well for 
a broad power-law distribution of field photons n'^ ex e'^"''' 
for different photon spectral indices a^, and readily explains 
the significant increase in the mean inelasticity of incident 
protons (or neutrons) from (Kp^) ~ 0.2 for steep photon 
spectra n' jj(e') with a^ > 1, to (Kp^) -^ 0.6 for hard spec- 
tra wifli a^, < 1 (Miicke et al. 1999). 

Calculations of the spectra of secondary particles and pho- 
tons (t^, 7, e) are done in the (5-function approximation, 
assuming that the probabilities for producing pions of dif- 
ferent charges (tt^, 7r+ and tt^) are equal for the multi-pion 
interaction channel. To correctly apply the (5-function ap- 
proximation, one has to properly take into account the differ- 
ent inelasticities of the multi-pion and single-pion production 
channels, which could explain the statement of Miicke et al. 
(1999) that the i5-function approximation does not work well 
in the case of hard photon spectra. 

2.5 Photon fields in the gamma-ray production region 

We consider 3C 279 as a prototype FSRQ. The broad-band 
radiation spectra produced by the relativistic ejecta in the jet 
of 3C 279 is generally thought to consist of a low energy syn- 
chrotron component and high energy Compton component. 
In 3C 279, the synchrotron component dominates at observed 
frequencies i' ^ 10^^ Hz. We assume that the v ^ 10^^ Hz 
synchrotron emission is cospatially located with the gamma- 
ray production region. This assumption is supported by ob- 
servations of correlated variability between the synchrotron 



and Compton components in some BL Lac objects and FS- 
RQs. The high-energy Compton emission consists of an SSC 
component and an external Compton (EC) component due to 
target photons that originate from outside the jet; for exam- 
ple from the disk, torus, scattered-disk radiation, reflected 
synchrotron radiation, etc. 

For calculations of n' j^(e') from 3C 279, we approximate 
the differential photon flux $(e) observed during the power- 
ful flare of 1996 (Wehrle et al. 1998) in the form of a con- 
tinuous broken power-law function, with power-law photon 
indices ai ^ 1.5, a2 — 2.45, and 0^3 ^ 1.6 at frequen- 
cies 1^ < 1^1 ^ 10^3 Hz, 1^1 < ly < 1^2 = 10^*^ Hz, and 
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i^ > 1^2, respectively. The lyFi, energy flux /(e) 
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of the synchrotron radiation reaches a maximum value of 
1.66 X 10^^"crgcm^^s~^ at 1/ = vi. The maximum flux 
of the Compton component is ~ 10 times larger than the 
peak synchrotron flux during the flare, and is reached at ~ 
500 McV. We suppose that the SSC component makes a sub- 
stantial contribution to the hard X-ray emission in 3C 279 
(Hartman et al. 2001), and therefore that a = 10. 

The redshift of 3C 279 is z = 0.538, and its luminos- 
ity distance is di^{z) = 1.05 x 10^^ cm for an f7,„ == 0.7, 
f2A ~ 0.3 cosmology with a Hubble constant of 65 km s^^ 
Mpc^^. The measured > 10^^ Hz flux during the 1996 flare 
implies that 5 > 4.9[£;p/i(GeV)°"^/[t„ar(d)]°-^^ from the 
77 transparency arguments. Thus (5 w 5 corresponds to the 
minimum Doppler-factor required for relativistic blobs in 3C 
279 to become transparent to absorption of the observed GeV 
photons in the radiation field inside the blob. 

Expressing r^ via the observed variability timescale ivar ~ 
1-2 days as observed by EGRET (Wehrle et al. 1998), one 
finds from equation (1) that 
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Here /_io = /(e)/(10-i"erg cm^^ g-i), and ^5 = (5/5. 
From the observed flux of synchrotron radiation from 3C 279 
with photon number index a = 1.5 at i^ <C 10^'' Hz, we de- 
duce B{G) ~ 20 r7"2/^[tvar(d)]"^/^ (Jj^"^^ from equiparti- 
tion arguments. The estimation of B from equation (2) gives 
B(G) ^ 32 t-^^{d)S^^, using L, = 3 x 10''^ ergs s-\ 
Lc/Ls = 10, and a = 10. 

The external radiation energy density in the blob rest frame 
is, using the expression for B derived from equipartition ar- 
guments, given by 

<,, ^ 20 (^) rj-^/'[tMd)r'^'sf'/' erg cm-^ . (5) 

For the calculations we assume rj = 1 and take a = 10. 



3 Calculations and Results 

We assume that the spectrum of the external UV radiation 
field arises from a Shakura/Sunyaev (1973) optically-thick 
accretion disk model that is scattered by broad-line region 
clouds. This disk model has flux density F{e) oc e^/"^ up to a 
maximum photon energy e„iaa; determined by the innermost 
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Fig. 1. Photomeson interaction energy loss time scale of protons, 
calculated for spectral fluxes observed from 3C 279 (see text) and 
tvar = 1 d, assuming 3 different Doppler-factors for the jet: 5 = 7 
(solid curve), 5 = 10 (dashed curve), and (5 = 15 (dot-dashed 
curve). The dotted and triple-dot-dashed curves are calculated for 
5 = 7 and 5 = 10, respectively, when p7 interactions with the 
synchrotron radiation field alone are considered. 



radius of the blackbody disk and properties of the central en- 
gine. We take rrieC^emax — 35 eV (Dermer and Schlickeiser 
1993). The energy density of the external radiation in the 
comoving plasma blob frame can be determined by noting 
that /(e }t 1)/ f{e <, 1) corresponds to the ratio Lec/^s of 
the directional external Compton and synchrotron powers as 
measured along the jet axis. 

Fig. 1 shows energy loss timescales (in the observer frame) 
of protons due to photopion production in a jet of 3C 279 cal- 
culated for 3 different Doppler-factors: 6 — 7 (solid curve), 
6 = 10 (dashed curve), and S — 15 (dot-dashed curve). 
The dotted and triple-dot-dashed curves show the photopion 
timescales corresponding to interactions with the synchrotron 
radiation only, for the cases S — 7 and (5—10. 

Fig. 2 shows the energy distributions of relativistic pro- 
tons A^p in the ejecta/blobs (for the same cases as in Fig. 1), 
which are calculated assuming power-law injection of rel- 
ativistic protons with number index ap = 2 on observed 
timescales At = 2 weeks. The total injection power of the 
protons Lp — 10''*(5~^ ergs s^^. Note that 10^* ergs s^^ is 
the characteristic apparent luminosity of the flares detected 
by EGRET from 3C 279. In calculations of iVp we take into 
account the photohadron interaction energy losses, as well as 
the escape losses of the protons in the Bohm diffusion limit. 
This limit is more restrictive than the condition that the gyro- 
radius of the highest energies of accelerated protons should 
not exceed the blob size (Hillas 1984). 

Fig. 3 shows the expected energy fluxes of the neutrinos 
produced by protons in Fig. 2. For the fluxes in Fig. 3, the 
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Fig. 2. The energy distribution of relativistic protons A^p calculated 
for the same jet Lorentz-factors 5 = 7, 10, and 15 and assumptions 
for the external and internal radiation fields as made in Fig. 1, as- 
suming a power-law injection of relativistic protons with number 
index Qp = 2 during At = 2 weeks with a power Lp = 10** J^* 
ergs s^"*". 



total number of neutrinos that could be detected by a 1 km^ 
detector during 2 weeks of the flare, using neutrino detection 
efficiencies given by (Gaisser et al., 1995) is 0.19, 0.12 and 
0.05 for the cases of 5 = 7, 10 and 15, respectively. It is 
important to note that if the external UV field is neglected, 
these numbers would be significantly reduced. We calculate 
the number of neutrinos that would be detected by a km^ 
array ranging from 1.7 x 10~^ to 3 x lO^'' for the fluxes 
in Fig. 3, which would not leave a realistic prospect for the 
detection of at least 2-3 neutrinos, which is required for a 
positive detection of a source. BL Lac objects, which have 
weak broad line regions and, therefore, a weak scattered ex- 
ternal radiation field, should consequently be much weaker 
neutrino sources. 

Because the properties of the external radiation field are 
insensitive to the value of the blob Lorentz factor, neutrino 
production in the quiescent state can be accurately estimated 
on the basis of properties of the external radiation field de- 
rived during the flaring state. For a typical duty factor of the 
flares ^1-2 months per year, and considering the additional 
neutrino production during the quiescent phase, we can rea- 
sonably expect that IceCube or other km^ array may detect 
several neutrinos per year from 3C 279-type blazar jets with 
(5^10. Allowing for the possibility that the nonthermal pro- 
tons may have an overall power larger than that for primary 
electrons improves the prospects for neutrino detection. 

It is important to notice that the external photon field re- 
duces very significantly the characteristic energy of protons 
for effective photopion production, compared with the case 
of the internal synchrotron radiation field. In fact, acceler- 
ation of protons only to E-p ^ 10^'' — 10^^ cV is required 
for efficient neutrino production through photomeson inter- 
actions, as can be seen from Fig.l. This effect permits us 
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Fig. 3. The fluxes of neutrinos expected due to photomeson inter- 
actions of protons shown in Fig. 2 with the external (UV) and in- 
ternal (synchrotron) photons in the jet of 3C 279 for Dopper factors 
5 = 7 (solid curve) 5 = 10 (dashed curve), and 5 = 15 (dot-dashed 
curve). The dotted and triple-dot-dashed curves show the calculated 
neutrino fluxes if the external photon field in the jet is neglected. 



to predict that a km^ array will detect high-energy neutri- 
nos from FSRQs without requiring acceleration of protons to 
ultra-high energies Ep ^ 10^* cV. 
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